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Abstract 

Heavy quark production is studied at the high energy linear e + e~ collider 
(LC) TESLA both in its nominal and Photon Collider (PC) mode. Leading 
order cross-sections are calculated for the production of heavy quarks, e + e~ — > 
e + e~ Q(Q)X, at high transverse momenta. The sensitivity of this process to the 
gluon content in the photon is studied. 



An e + e~ linear collider with a centre of mass system (CMS) energy in the range of 
500 to 1000 GeV is the prime candidate for the next high energy accelerator project 
after the LHC Recently TESLA [p]] has completed a Technical Design Report for such a 
collider. For the study of a future electron-positron Linear Collider (LC) it is important 
to examine the physics potential of its main and possible additional options. Such an 
option is the Photon (or Compton) Collider (PC) in which high energy real photons can 
be obtained by backscattering photons from a laser beam on the electron or positron 
beam 0, [|. This way an excellent tool for the study of 77 collisions at high energies 
can be provided. 

In high energy e + e~ collisions the hadronic final state is predominantly produced 
in 7*7* interactions where the virtual photons (7*) are almost on mass shell. Their 
scattering can be effectively described with the interaction of real photons with energy 
spectra given by the Weizsacker- Williams (WW) approximation [[|. Such a spectrum 
is generally quite soft, as shown in Fig. 1. In contrast, a Photon Collider, based on 
Compton scattering, provides beams of real photons which are much harder than those 
from a WW distribution (see Fig. 1). Furthermore these photons can be produced 
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in a definite polarization state. While the energy spectrum of the PC at the high 
energy end can be well calculated based on a few parameters of the laser and incident 
electron beam, the lower energy end depends much more on the details of machine and 
operation. Therefore we limit ourself to study only the high energy peak of 0.6 < y = 
E 1 /E e < 0.83 ||. The spectra used in this paper are described in the Appendix. 

The main goal of this study is to compare the LC and PC potential for probing 
the gluon distribution in the photon via heavy quark production. This analysis does 
not make use of polarization. The measurement of the process e + e~ — > e + e~Q(Q)X is 
an important test of QCD. It is well known that this process is promising for probing 
the structure of the photon, see e.g. || [7j and ||, where a related topic is considered. 
Heavy quarks can be produced in 77 collisions through three mechanisms. Direct 
(DD) production occurs when both photons couple directly to a QQ pair. In single 
resolved photoproduction processes (DR) one of photons interacts via its partons with 
the second photon. When both photons split into a flux of quarks and gluons, the 
process is labelled a double resolved photon (RR) process. 

There are various schemes proposed for calculating cross-sections for processes in- 
volving heavy quarks, for the latest overview see ref . Two standard schemes will be 
used in this analysis Q. They differ by the number of quark flavours which correspond 
to partons of the photon (active quarks). In the case of the massive quark scenario, the 
so called Fixed Flavour Number Scheme (FFNS), the photon "consists" only of light 
quarks and gluons, which may interact, and massive heavy quarks can be only created 
in a hard process e.g. via gluon-gluon fusion. The massless quark scheme, Zero-mass 
Variable Flavour Number Scheme (ZVFNS), considers apart from the gluons and u, d 
and s quarks also heavy quarks as active quarks, which are all treated massless. All the 
partonic reactions in both schemes contributing in the LO calculation to heavy quark 
production are shown in table 1. Present available data on charm quark production in 
77 processes is still not precise enough to allow to discriminate between the predictions 
of these two production schemes. Besides, these theoretical approaches are expected 
to be in turn appropriate for describing different ranges of the hard scale (/i) of the 
process: The massive and massless schemes should be applied below (/i < mq) and 
above (/i ^> tuq) the mass scale, respectively. For the charm quark and for transverse 
momentum of the partons larger than e.g. 10 GeV, we expect to be closer to the 
region of applicability of the ZVFNS scheme than the FFNS scheme. For completeness 
we performed the calculation in both schemes which can be treated as an estimate of 
the theoretical uncertainty. 

We calculate the LO QCD cross-sections for c and b quarks produced with large 
p T at the e + e~ LC colliders at energies of 300, 500 and 800 GeV, respectively, and at 
the PC based on the corresponding (e + e~) LC with these energies. For comparison we 
study also the corresponding cross-sections for the LEP collider at 180 GeV. To test 
the sensitivity of the considered processes to the gluonic content of the photon we use 
two different parton parametrizations for the real photon: GRV 1TT1 and SaSld [IIT 



x a recent attempt to construct a new composite scheme can be found in 
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We choose the GRV('92) parametrization because in its construction heavy quarks are 
treated as massless above the heavy quark threshold region (W 3> 2m<g), and moreover 
recent experimental results IT^j indicate that GRV gluon distribution is close to the 
data. The choice of the SaSld parametrization is motivated by the fact that it provides 
results in agreement with expectations of Gribov factorization and the QCD sum rule 



for the structure of the photon [13]. Both these parametrizations are based on QCD 



fits to photon structure function data measured in cy collisions at e + e experiments 



14, 15 1 . They have different assumptions for the gluon content, which is only weakly 
constrained by the data. The GRV and SASld distributions both start the evolution 
from a small scale, /ijj = 0.25 GeV 2 and 0.36 GeV 2 , respectively, and consequently 
both parton densities predict a rise of the gluon density at small x. Different treatment 
of valence quark distributions in the vector mesons leads to a larger gluon component 
at small-x 7 of the photon for the GRV set compared to the SASld one. 

As was mentioned above, in the massive (FFNS) calculations the number of active 
flavours (Nf) is taken to be 3. In the massless (ZVFNS) scheme it varies from 3 to 5 
depending on the value of the hard (factorization, renormalization) scale \i. The c (b) 
quarks are included in the computation provided that /i > m c (nib) with m c = 1.5 GeV 
(nit, = 4.5 GeV). When charm quark production is calculated in the ZVFNS scheme 
the bottom quarks are always excluded, hence Nf = 3 or Nf = 4. Also the QCD energy 
scale Aqcd, which appears in the one loop formula for the strong coupling constant 
a s , depends on the number of active flavours. We take this scale, denoted as Aq CD) to 
be: 

A* CD = 232 A* CD = 200 A 5 QCD = 153 MeV (1) 



as m 



TUfl. The hard scale in the calculation of the cross-section \i is taken to be the 



transverse mass of the produced heavy quark my = \JmQ + p\. 

In the calculation we include direct (DD), single (DR) and double resolved (RR) 

photon processes. We study , with r\ =| In the rapidity of the produced 

heavy quark and px its transverse momentum. The results for the charm quark pro- 
duction presented in Figs. 2 and 3 were obtained in the ZVFNS scheme for both 
types of initial photon spectra. The calculation was performed using the GRV parton 
parametrization. In Fig. 2 a fixed energy a/s=500 GeV and a, pr =10 GeV for the 
charm quark was assumed. Both resolved photon contributions, DR and RR, to the 
process e + e~ — > e + e~Q(Q)X are dominated by reactions initiated by gluons, especially 
for a PC at high rj values. The dominance of the gluon initiated over the quark initi- 
ated processes is larger for the PC spectrum; it is also larger in the FFNS scheme (not 
shown) compared to the ZVFNS one. 

Fig. 3 shows a comparison of the direct (DD), single resolved (DR), double resolved 
(RR) contributions and their sum for WW and LASER spectra for various CMS en- 
ergies of the e + e~ colliders, for charm quarks with a pt of 10 GeV. An interesting 
pattern is observed. In the case of an e + e~ LC with a WW photon spectrum either 
the direct production of c dominates, or resolved and direct photon contributions are 



3 



found to be of the same size. The DR and RR contributions increase with increasing 
CMS energy of LC. Nevertheless in the considered range of the LC energies, they do 
not play a dominant role for heavy quark [] production in the e + e~ mode. The contrary 
is observed for a PC: charm quark production is always dominated by resolved photon 
interactions. The direct contribution is limited to high absolute values of the rapidity 
distribution. For a pr = 10 GeV this rapidity region is outside the range of \rj\ < 2 
which could be covered by a typical detector. Lowering the px of the observed charm 
quarks amplifies the contribution of the resolved process to the total cross-section. The 
direct contribution becomes even less important with increasing CMS energy. Thus the 
charm production cross-section is much more sensitive to the parton distribution of the 
photon for a PC compared to a LC. Since for the PC option the resolved photon contri- 
butions are dominated by the processes involving gluons (Fig. 2) this offers an excellent 
tool for measuring the gluonic content of the photon. We note here also that an ex- 
tention of the experimental reach to observe charm quarks at rapidities beyond \rj\ — 2 
would allow to kinematic range where the gluon contribution dominates most 

strongly over other resolved photon contributions. 

For the ZVFNS scheme and the PC option we find that the processes initiated by 
c quarks constitute about 52-67% of the total cross-sections presented in Fig. 3, and 
about 99% of the contribution which arises due to the processes initiated by quarks 
(not shown). Since the heavy quark content of the photon is closely related to the gluon 
density in the photon, the heavy quark production discussed in the paper exhibits also 
sensitivity to the gluonic density via the heavy quark content of the photon Ul6"f . Since 



for the LC option the cross-section has a large contribution from the direct process, this 
effect should be most visible in data from a PC, which are dominated by the resolved 
photon processes instead. 

An important feature of the results for cross-section seen in Fig. 3 is the increase of 
the resolved photon process contributions, and consequently also an growing sensitivity 
to gluons, with increasing CMS energy. This arises from the fact that at higher energies 
we probe regions of small Bjorken-x 7 values. The minimal x 7 value that could be 
reached in LC for p T = 10 GeV varies from ~ 0.004 for y/7s = 300 GeV to ~ 0.0006 
for -y/s = 800 GeV (~ 0.01 in LEP). The detector cut \rj\ < 2 applied in this analysis 
increases those values to ~ 0.006 and ~ 0.0018, respectively. An extra cut on the 
LASER photon spectrum (0.6 < y < 0.83) in the PC mode shrinks the x 7 accesible 
range even further. The numbers read then: ~ 0.008 for = 300 GeV and ~ 0.0023 
for t/s = 800 GeV. Therefore the RR contribution grows weaker with the CMS energy 
in the PC case comparing to the LC case. For the DR case in the PC mode the 
0.6 < y < 0.83 cut puts an additional, kinematic constraint. Since one of the photons 
interacts directly with E y > 0.QE e the energy of the produced heavy quarks must grow 
with E e . If we choose them to have relatively small pt equal to 10 GeV it is more 
difficult to produce them in the central rj region where E ~ px- Hence instead of the 
increase with the energy we see a slight decrease of the DR contribution around rj = 0. 
Due to the same kinematic effect DD events cannot be produced close to r\ — 0. We 

2 A similar effect is observed for b quark production (not shown) 
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note as above that reaching in experiment values of rapidity higher than \r]\ > 2 would 
be of high importance. We then could explore kinematic region of strongest resolved 
contributions to the cross-sections. 



The results discussed so far are not significantly affected by the choice of scheme for 
the heavy quark calculation as is shown in Fig. 4 for the c-quark rapidity distribution: 
for both the FFNS and ZVFNS scheme the cross-section is much larger for the PC 
than for the LC case. 

More details of the comparison of the FFNS and ZVFNS scheme results can be found 
in Figs. 5-9. The importance of the low x 7 contributions is demonstrated in Fig. 5, 
presenting the 4^- distributions of the charm production at a PC with J~s^ e = 500 



dx-y 

GeV (for the GRV parton parametrization), on a linear and logarithmic rc 7 scale. 
The cross-section has been integrated over the region p T > 10 GeV and \rj\ < 2, 
realistic for a possible detector at a PC. The direct contribution is included into the 
bin with the highest rr 7 . The double resolved part is counted for twice, once for each x 7 
value. The DR contributions are peaked at small x 7 , while the RR ones have a flatter 
distribution. The expected differences between both schemes used in the calculation 
are visible: the ZVFNS scheme gives a larger contribution of DR and RR processes, 
especially for a; 7 >0.2. The maximum in the DR distribution occurs at x 7 close to 
0.01, independently on the scheme. The RR process have a broad distribution between 
x 7 ~ 0.03 — 0.05 and 1.0, and the ZVFNS scheme result is larger by more than one 
order of magnitude compared to the FFNS one in this whole region. 

The sensitivity of the charm quark production process to the gluon distribution 
is studied further by comparing the predictions obtained using two different parton 
parametrizations for the photon. In Fig. 6 (7) the ratio of the relative difference 

of the cross-sections is presented, obtained using the GRV and SaSld parton 

parametrizations in the ZVFNS (FFNS) scheme. As expected the cross-section ob- 
tained using the PC photon spectrum (LASER) leads to a larger sensitivity than the 
WW spectrum for a given energy of the e + e~ collider: the difference between results of 
these two parametrizations shown is 5-20% for a WW and 30-50% for LASER photon 
spectrum. Hence a PC collider has a larger sensitivity for testing the photon structure 
and especially its gluon density. It is worth mention that these results indicate however 
that an initial study of the partonic structure of real photon can be already made with 
an e + e~ LC itself. 

Next we show some detailed characteristics of charm quark production at a PC. 
Fig. 8 and 9 show the relative difference of x 7 distribution calculated with the GRV 

and SaSld parametrizations, for the J^— cross-sections introduced before, for \rj\ < 2. 

Results for p T > 10 GeV are shown in Fig. 8 (for both linear (A) and logarithmic (B) 
scale) and for p T > 5 GeV in Fig. 9 (linear scale only). 

Table 2 shows in detail how the gluon initiated resolved processes dominate the 
charm quark production. For the case of y^s = 800 GeV they contribute 76% and 
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of the integrated cross-section calculated in the FFNS and ZVFNS scheme respectively. 
Further details on contributions of the gluon and charm quark initiated processes can 
be found in Table 3. 

Up to now we presented only the results for c quark production. The corresponding 
b quark production has all the features listed above though the difference between a PC 
and a LC was found to be smaller. All calculated cross-sections for bottom production 
are found to be even more sensitive to the shape of gluonic content of the photon than 
the corresponding ones for charm quark production, but the cross-sections are smaller. 

Next we calculate the number of events with heavy quarks, c and b, assuming an 
e + e~ integrated luminosity of 300 fb" 1 , which could be achieved at a high luminosity 
LC with one year of running. The corresponding luminosity of the PC for the region 
0.6 < y < 0.83 is about a factor 3 to 4 lower [17]]. In the FFNS scheme the event 
rates for charm quark production at a LC for the different CMS energies 300-800 GeV 
amount to 1 — 4 ■ 10 6 events for p? > 10 GeV; the corresponding numbers at a PC 
are 9 — 25 • 10 6 events. Hence these calculations predict a high number of the heavy 
quarks, c and b, produced at the considered CMS energies of the e + e~ collider of 300, 
500 and 800 GeV for both the LC and PC options. 

In practice, charm quarks with a > 5 GeV produced e.g. via D* decays can be 
detected in a generic LC detector without dedicated detectors in the rapidity range 
of \r}\ < 1.5 — 2. This restricts the lowest reachable x 7 values to 0.0011 and 0.0023 
for a minimum px of 5 and 10 GeV, respectively, for a PC collider at 800 GeV. The 
cross-sections for c and b production in this phase space region are presented in Table 
4. The corresponding event rates (for FFNS scheme leading to smaller rates) are given 
in Table 5. It shows that the event rate at a PC is still larger than the one at a LC, 
but the difference is reduced to a factor 2 to 4. Hence some of the advantage of the 
PC over the LC is lost due to the charm quark production at large rj in case of the PC 
(see Fig. 4), which will go undetected with the presently planned detectors. 

To measure the charm content of the photon e.g. D* mesons can be selected to 
reconstruct the kinematics of the scattering process. The efficiency including fragmen- 
tation fraction and branching ratios, is typically around a few times 10 -3 |T%fl . Hence 
the number of reconstructed events per year with charm quark with px > 10 GeV will 
be approximately a few thousand for the LC and about ten thousand events for the 
PC. The numbers are about 5 times larger for a reduced cut of pt > 5 GeV. These data 
will clearly allow for precision measurements of charm quark production and further 
on the gluon distribution in the photon. The statistical precision of the measurements 

of dp't'dr) w ^ n e -S- ten bins will be approximately 5% at the LC and a few % at the 
PC. The event rates in the ZVFNS scheme can be larger by a factor 3, compared to 
the FFNS scheme. 

In conclusion, the calculated cross-sections for heavy quarks (c and b) production 
in two photon collisions show a higher sensitivity to the parton distributions of the 
photon in case of a Photon Collider compared to a e + e~ Linear Collider. This does 
not depend on the particular scheme used to calculate the heavy quark cross-sections. 
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Since the resolved photon contribution is to a large extent dominated by gluon induced 
processes, especially for a high energy PC, we conclude that heavy quark production 
provides indeed a sensitive probe of the gluon content of the photon. Combining the 
above features with the larger cross-sections achieved at energies of the e + e~ collisions 
at a PC favours this option for future photon structure research. A high luminosity 
e + e~ collider which drives the PC collider will however be essential. 



Appendix 

The following formula of the Equivalent Photon Approximation is used: 
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In case of the PC mode we use the original energy spectrum of unpolarized photons 
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where k is a parameter giving the restriction of y value: y< jt^- It is argued that 
the optimal value of k is 4.83, which gives a cut of y, y ma x = 0.83. We have chosen 
these values for this analysis. Because the part of the spectrum with y < 0.6 is very 
sensitive to the technical parameters of the PC such as the size of the beam we used 
only its high energy peak of 0.6 < y < 0.83 B. 
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DD 


DR 


RR 


FFNS 


7 + 7 ^Q + Q 


g + l^Q + Q 


g + g^Q + Q 

q + q-^Q + Q 


ZVFNS 




g + l^Q + Q 
Q(Q)+i^Q(Q) + g 


g + g^Q + Q 

q + q^Q + Q 
Q+Q^Q+Q 
Q + Q(Q + Q) -> Q + Q(Q + Q) 
Q(Q)+q(q)^Q(Q)+q(q) 

Q(Q) + g^Q(Q)+g 



Table 1: Parton reactions contributing to heavy quark production in 77 collision in 
the massive (FFNS) and massless (ZVFNS) schemes. 
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Figure 1: The photon spectra f 7 (y) as a function of the energy fraction of the photon. 
A comparison is shown of the spectrum of photons radiated by the initial electron 
according to the Weizsacker- Williams (WW) approximation with the spectrum of (un- 
polarized) photons produced in the backscattering Compton process on the laser beam 
(LASER) for k = 4.83, see Appendix for definitions. For the LASER case the high 
energy peak y > 0.6 is indicated. 
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Figure 2: Comparison of parton processes which involve gluons with other partonic con- 

tributions to the resolved (DR+RR) part of the cross-section ^-g ^ e + e — > e + e c/ cX 

for y/s = 500 GeV. Results are obtained for charm quarks with transverse momentum 
Pt = 10 GeV as a function of rapidity rj (in the e + e~ CM frame) in the ZVFNS scheme 
with the GRV parton parametrization using: A) the WW spectrum, B) the LASER 
spectrum. 
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Figure 3: The cross-sections -j-j 



dpfr^drj 



(e + e — > e + e c/cX) at various LC energies as a 



function of the c/c rapidity 77 in the e + e~ CM frame at pr = 10 GeV presented for two 
types of photon spectra (WW and LASER). Separately the DD, DR, RR and their 
sum are shown. For a comparison also estimates for LEP 2 (WW, ^=180 GeV) are 
given. The GRV parton parametrization and the ZVFNS scheme are used. 
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Figure 4: The cross-section 



dy 



e + e c/cX) at ^/s = 500 GeV as a 



function of i] for = 10 GeV. In the calculation the GRV parton parametrization is 
used. A comparison is shown between ZVFNS and FFNS schemes, and between the 
WW and LASER photon spectra. 
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Figure 5: The cross-section J^— (e + e — > e + e c/cX) as a function of x 7 for yfs = 500 

GeV, integrated over p T > 10 GeV and < 2 for two types of schemes, ZVFNS and 
FFNS. The LASER spectrum and the GRV parton parametrization are used; A) linear, 
B) logarithmic x 7 scale. 
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function of 77 for various LC energies and for px = 10 GeV. Results are obtained in the 
ZVFNS scheme for two photon spectra: A) WW, B) LASER. 
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Figure 7: As in Fig. 6, for the FFNS scheme. 



13 



DR 



RR 



DD+DR+RR 



0.6 



- 1 1 1 1 1 1 1 

: ZVFNS - 


i 1 : 


'- FFNS - 

: _ a") l "l 




-.i/i.i, 


1 1 1 1 



0.6 


: ' 1 ' 1 


I'M; 


0.6 


: ' 1 ' 1 


■ 1 ' 1 ' . 


0.4 


~_ 


1 


0.4 


j- i -J L 




0.2 






0.2 




- 
















0.2 


'■ , i , i 


1 , 1 , : 


-0.2 


: , I , I 


. 1 , 1 , : 



0.2 0.4 0.6 0.8 

DD+DR+RR 
_ | — i 1 1 1 Miij 




da 



(e + e — > e + e c/cX) as a 



Figure 8: The ratio ^ RV Q^y^ ld °f t ne cross-section 

function of x 7 at y/s = 500 GeV for > 10 GeV and \rj\ < 2. A comparison is shown 
of the two schemes for a LASER spectrum on A) linear and B) logarithmic x 7 scale. 
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Figure 9: As in Fig. 8, for > 5 GeV and linear scale only. 
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^ = 300 GeV 


3.1 


40 


5.6 


0.34 


5.5 
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17 
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1.1 


43 
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0.37 
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0.43 


42 


3.8 


0.41 


7.5 


7.9 


24 
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Table 2: Cross-sections (in p6) for charm quark production with p T > 10 GeV and 
\rj\ < 2 (0.6 < y < 0.83), calculated in the massless ZVFNS scheme with the GRV 
parton parametrization. 



p T > 10 GeV 


ZVFNS 


FFNS 




cr (f)/ (7 


a(c)/a 


a(g + c)/a 


°(9)/v 


yfe = 300 GeV 


0.70 


0.58 


0.97 


0.92 


yfe = 500 GeV 


0.74 


0.61 


0.99 


0.97 


yfe = 800 GeV 


0.76 


0.66 


0.99 


0.98 



Table 3: An estimate of the partonic contributions to the charm quark production with 
p T > 10 GeV and \rj\ < 2 (0.6 < y < 0.83), calculated in both the massless ZVFNS 
and massive FFNS schemes with the GRV parton parametrization. The sum of the 
contributions which involve gluons (a(g)), charm quarks (cr(c)) or both (a(g + c)) in 
the initial state are compared to the total cross-section a. 
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p T > 10 GeV 


p T > 5 GeV 




c/c 


b/b 


c/c 




LC 


PC 


LC 


PC 


LC 


PC 


= 300 GeV 


4.0 (5.4) 


42 (110) 


0.27 (0.22) 


8.4 (23) 


25 (39) 


170 (490) 


yfe = 500 GeV 


6.0 (8.6) 


45 (120) 


0.46 (0.94) 


10 (32) 


35 (58) 


180 (600) 


yfe = 800 GeV 


8.2 (13) 


47 (150) 


0.73 (1.6) 


13 (47) 


47 (83) 


200 (770) 



Table 4: Cross-sections (in pb) for heavy quark production with \q\ < 2, at a LC and a 
PC (0.6 < y < 0.83), calculated in the massive FFNS (ZVFNS) scheme with the GRV 
parton parametrization. 





p T > 10 GeV 


p T > 5 GeV 




c/c 


b/b 


c/c 




LC 


PC 


LC 


PC 


LC 


PC 


yfe = 300 GeV 


1.2 


3.9 


0.08 


0.80 


7.5 


16 


^ = 500 GeV 


1.8 


4.1 


0.14 


0.97 


11 


17 


yfe = 800 GeV 


2.5 


4.4 


0.22 


1.2 


14 


19 



Table 5: Event rates (in units of 10 6 ) for heavy quark production with \r)\ < 2 at a 
LC and a PC (0.6 < y < 0.83), calculated in the massive FFNS scheme with the GRV 
parton parametrization. The e + e~ LC luminosity was assumed to be 300 fb -1 . 
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